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Current-carrying and superconducting systems can be treated within density-functional theory if 
suitable additional density variables (the current density and the superconducting order parameter, 
respectively) are included in the density-functional formalism. Here we show that the corresponding 
conjugate potentials (vector and pair potentials, respectively) are not uniquely determined by the 
densities. The Hohenberg-Kohn theorem of these generalized density-functional theories is thus 
weaker than the original one. We give explicit examples and explore some consequences. 
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Today, density-functional theory (DFT) fll IS an in- 
dispensable tool for the investigation of the electronic 
structure of matter in atomic, molecular, or extended 
systems. The theory rests on the celebrated Hohenberg- 
Kohn (HK) theorem ||], which guarantees that the (v- 
representable) ground-state density n(r) uniquely de- 
termines the ground-state many-body wave function 
'00 (ri, . . . , r^v). This theorem on its own is a very pow- 
erful result, but in the original formulation of DFT [^|J|] 
one can prove even more: the external potential t;(r) 
(e.g., the nuclear charge distribution in a molecule or a 
solid), too, is a functional of the density, and is unique 
up to an additive constant. Since this external poten- 
tial in turn determines all eigenstates of the many-body 
Hamiltonian, this implies that all observables (and not 
only ground-state ones) are functionals of the ground- 
state density. 

Following original ideas of von Barth and Hcdin Q , it 
has recently been shown by Eschrig and Pickett [p| and 
by the present authors §] that in spin-DFT (SDFT) the 
situation is not that simple: while the wave function is 
still uniquely determined by the spin densities rif (r) and 
nx(r), the external potentials i'-f(r) and ^^(r) [or w(r) 
and B(r)] are not. This implies that SDFT function- 
als are not always differentiable, and has far-reaching 
consequences for the construction of better exchange- 
correlation (xc) functionals, and for applications to sys- 
tems such as half- metallic ferromagnets [^,^ . 

SDFT is not the only instance at which the original 
HK theorem has been generalized. In the present work 
we extend the analysis of Ref. |^ to two other gener- 
alizations of DFT, namely current-DFT (CDFT) [@,§. 
and DFT for superconductors |9|-[T^. The discovery of 
nonuniqueness in these generalized DFTs deepens our un- 
derstanding of the respective xc functionals and flags a 
warning signal to all-to-immediate generalizations of the 
original HK theorem to more complex situations. 

The basic physics of nonuniqueness is simple. When 
a sufficiently small change in one of the external fields 



does not change the corresponding density distribution, 
the corresponding susceptibility vanishes. The search 
for, and the interpretation of, nonuniqueness in DFT is 
thus guided by investigations of the circumstances under 
which some response function becomes zero. 

We first consider current-carrying systems. The ap- 
propriate formulation of (nonrelativistic) DFT is CDFT 
which is based on the many-body Hamiltonian (in 
atomic units, i.e., fi = e = m = 1) 

H^f + U + J d^rh{r)[v{r) - + ^ J (iVjp(r)A(r) 

+ ^ j d\fi{Y)A^{r) + y"dVm(r)B(r), (1) 

where B(r) = V x A(r) is the magnetic field, ^(r) the 
electrostatic one, and T and U denote the operators for 
kinetic energy and particle-particle interaction, respec- 
tively. 

The basic variables of CDFT, in terms of which the 
entire ground-state physics of the current-carrying many- 
body system is described, are n(r), m(r), and jp(r), the 
ground-state expectation values of the particle density 
operator n(r) = ^],.(r)\E'cr(r), spin magnetization op- 
erator m(r) = (l/2c) /3 ^a(r)o'^^(r), and (param- 
agnetic) current density operator 

J^-W = - (V*Ur))vI/.(r)] , (2) 



where the ^^(r) are field operators and a is the vector 
of Pauli matrices. 

According to the CDFT version of the HK theo- 
rem these densities uniquely determine the ground-state 
many-body wave function. However, in striking contrast 
to conventional 'density-only' DFT they do not uniquely 
determine the potentials A(r), B(r), and w(r): it is pos- 
sible to find different vector and scalar potentials which 
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yield the same ground state, and consequently the same 
densities n(r), m(r), and jp(r). 

Before delving into a general characterization of such 
potentials, we present a simple example that clearly dis- 
plays the problem. We consider an atom subjected to 
a uniform magnetic field B — Bz, where z is the unit 
vector along the z-axis. Ignoring spin-orbit interactions, 
the Hamiltonian can be written as 



8c2 



E — 



2c 



(3) 



where ', 



±,i 



and L2 and Sz are the z-components 



of the orbital and spin angular momentum operators L 
and S. Z is the atomic number, specifying the exter- 
nal potential. Both and Sz are constants of motion, 
hence the (nondegenerate) ground state of H is also an 
eigenstate of Lz and 6*2, with eigenvalues and mg 
respectively. The ground-state energy is Eq. 

Consider now the same system of Z electrons being 
subjected to a different (but still uniform) magnetic field 
B' = B'z and the external potential 



Z 
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v\r) = ---—{B''-Byi. 
r 8c^ 

The Hamiltonian of this system is 

H' = H+ ^{Lz + 2Sz){B' - B). 



(4) 



(5) 



Thus, we immediately see that the ground state of B 
(or, for that matter, any simultaneous eigenstate of i?, 
Lz-, and Sz) is also an eigenstate of B' with eigenvalue 
E' = Eq^ (l/2c)(rnL -f 2ms){B' - B). Furthermore, if 
the difference B' — B is not too large, this eigenstate will 
be the ground state of B'; the qualitative condition for 
this to happen is that E' — Eo « Eg where Eq is the 
energy gap between the first excited state and the ground 
state of B. Thus, we have succeeded in constructing two 
different sets of potentials, A = (B x r)/2 and v, and 
A' = (B' X r)/2 and v' , that yield the same ground state. 

Let us now consider the question from a more general 
point of view. Let A' = A + AA and v' ^ v + Av be 
vector and scalar potentials that are supposed to yield the 
same ground state ipo as A and v. A necessary condition 
for this is that tAo satisfy the eigenvalue equation 



fiAv + ArfiAA^ + -jpAA 
2c'' c 



V^o^ASV-o, (6) 



where we neglected, for simplicity, the spin-degrees of 
freedom, because the nonuniqueness associated with 
them is already discussed in Refs. ||] and ||]. 

The general problem at hand is thus to find a linear 
combination of the density operators n(r) and jp(r) that 
has ijjo as eigenf unction. This problem is not easily solved 



in general. It is easy, however, to obtain a particular so- 
lution of Eq. (|^) if one can find a linear combination of 
the density operators that is a constant of motion. The 
ground state of B is automatically an eigenstate of such 
a constant of motion, and Eq. (^ is satisfied. By mak- 
ing the coefficients of the linear combination sufficiently 
small we can always ensure that ipo remains the ground 
state of the Hamiltonian with the new potentials (assum- 
ing of course that the spectrum of B has a gap between 
its ground state and first excited state) . This is the same 
prescription employed in Ref. to construct examples 
for nonuniqueness in SDFT. In the terminology of that 
reference nonuniqueness arising from such constants of 
motion is referred to as systematic nonuniqueness. 

As a trivial example of this procedure consider the 
constant of motion N = J d^rn(r). The existence of 
this constant of motion tells us that At;(r) = const, 
AA(r) = is a solution of Eq. (|). This is the well 
known nonuniqueness of the scalar potential with respect 
to the addition of a constant. Consider now the less triv- 
ial example 



dV(z X r) • jp(r). 



(7) 



which is a constant of motion in any system that is in- 
variant under rotations about the z-axis. Comparing this 
with Eq. (^) we immediately see that AA(r) = AB{z x 
r)/2 and A'i;(r) = -[AA(r)]2/(2c2) with AB = const is 
indeed a solution of the posed problem. This is, of course, 
nothing but a more formal derivation of the elementary 
example discussed above. 

Another way in which nonunicjucness can arise is by 
adding an operator to the Hamiltonian which, although 
not a constant of motion, happens to have eigenvalue 
zero on the ground state. This was called acciden- 
tal nonuniqueness in Ref. |^. To give an example 
in CDFT, let no-(r) and jpcr(r) denote the exact spin- 
resolved ground-state density and paramagnetic current 
of a two-electron system, such as the Be atom, in the 
presence of external vector and scalar potentials. For 
sufficiently small external fields these densities must arise 
from the single particle orbitals (pi{r) and <fi{r) that are 
the lowest energy solutions of the spin-dependent Kohn- 
Sham (KS) equations 

(8) 

where Ascr(r) and Us(T(r) are the KS potentials, defined, 
as usual, in terms of the external, Hartree, and exchange- 
correlation potentials. The relation between the densities 
and the single-particle orbitals is ^^(r) = \(pcr{r)\'^ and 



jp^(r) 7i<^(r)V(/)^(r), 



(9) 



where <^cr(r) is the phase of the complex orbital (pa-{j^) 

ill. 
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Eq. (@) can be rewritten in the form 



(10) 



which has the solution Aso-(r) — 
-cjp„{r)/n„{r) and 



-cV^„(r) 



(r) 



(r) 



(11) 



To determine whether these are the only potentials 
that reproduce the given densities na{r) and jpo-(r) we 
assume the existence of a second such set of potentials, 
A',^(r) = Asa{r)+AA,^{r) and u^^) = w^(r)+Ai;,.^(r). 



By substituting these back in Eq. ( 
real and the imaginary parts we o 



lOD and separating the 
Dtain 



(r) + ±AA^ 



(12) 



and V • [tIct (r)AAs(j(r)] — 0. This last equation follows 
more directly from the application of the continuity equa- 
tion to the real solution of (p^). Its general solution is 
AAs(j(r) = V X Qcr(r)/no-(r) where Qcr(r) is an arbitrary 
vector field. Hence, 



(r) 



1 
2? 



and 



A;,(r)=A,,(r) 



V X Q^(r) 
na{r) 



V X Q^(r) 



(13) 



(14) 



By construction, Aso.(r) and Uso-(r) are the potentials for 
which ^^(r) and jpCT(i") are grrownd-siate densities. If Qa 
is sufficiently small and if the KS system at Qo- = has 
an energy gap separating the first excited state from the 
ground state, ^a-(j) will remain the ground state in the 
potentials A',^(r) and v',^{r). Thus, Eqs. (|l|) and (|l|) 
provide a vivid and nontrivial example of nonuniqueness 
of the KS potentials of CDFT. 

Next, we turn to the superconducting case. Here the 
underlying many-body Hamiltonian is 



H = f + U + J dVn(r)[w(r) -/i] + J d^r m(r)B(r) 

- J d^'' J d^r'[xir,r')D* {r,r') + H.c], (15) 

where the expectation value of the pair operator 
x(r,r') = 'I'|(r)5'|(r') is the superconducting order pa- 
rameter and -D(r, r') the corresponding pair potential. 
The phonon-induced interaction term of Rcf. p| can be 
added to H without changing our conclusions. 



As above, we now assume that the densities n(r), m(r) 
and x(r, r') can also be reproduced in different fields v' = 
?; Ad, B' = B -I- AB and D' = D + AD. The equation 
obeyed by Av, AB and AD is 



y (fr nAv + mAB- J (fr'{xAD + H. 



c.) 



^0 



= A£;^o. (16) 

At this stage we already see a first nontrivial difference 
to the case of CDFT and SDFT: due to the presence of 
the pair operator x'm H the particle number operator N 
is not a constant of motion, and we are not free to add an 
arbitrary constant to the external potential v{r). In other 
words AD = 0, AB = and Av = const is not a solution 
of Eq. dH) for a given -00- DFT for superconductors 
thus does not suffer from the most basic nonuniqueness 
of all, that with respect to the additive constant in the 
electrostatic potential. 

However, DFT for superconductors is not free of 
nonuniqueness. For a singlet superconductor the spin 
susceptibility vanishes at zero temperature ||l^. In the 
light of our physical characterization of nonuniqueness at 
the beginning of this paper we would thus expect some 
associated nonuniqueness. Indeed, this is bourne out by 
more detailed analysis. If B = Bi is spatially uniform 
and sufficiently weak not to break Cooper pairs param- 
agnetically, then B' = B + AB, where AB is also weak, 
uniform and parallel to z, has the same ground state, 
because under these circumstances = J (Pr rhz is a 
conserved quantity, i.e., the superconductor remains in a 
singlet state, with all electrons paired up. Consequently, 
the set of potentials {v, B, D} is not uniquely determined 
by the conjugate densities {n, m,x}. Since it is associ- 
ated with the constant of motion Mz , this is systematic 
nonuniqueness in the above sense. 

With these examples we end our list of explicit oc- 
curences of nonuniqueness in generalized DFTs, and now 
turn to a discussion of broader aspects of our findings. 

In early papers on both CDFT and DFT for su- 
perconductors |l^ one finds the statement that the cho- 
sen densities uniquely determine the corresponding po- 
tentials. As we have shown here, these statements are 
not accurate, and all that is determined uniquely is the 
ground-state wave function. Concerning consequences 
of this finding we refer the reader to the discussion we 
have given earlier of consequences of nonuniqueness in 
SDFT 1^. That discussion carries over almost literally 
to the case of current-carrying and superconducting sys- 
tems. However, we wish to stress particularly that for 
most applications of any DFT, including CDFT and DFT 
for superconductors, uniqueness of the ground-state wave 
function is sufficient, since no explicit use of the density- 
potential relation is made. A notable exception within 
CDFT is the recent work by Handy and Lee |jl^ , in which 
it is attempted to systematically construct exact CDFT 
potentials from given densities. This construction must 
be reexamined in view of our finding that the CDFT 
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potentials are not uniquely determined by the densities. 
Further exceptions are listed in Ref. 1^. 

Another important consequence of nonuniqueness in 
DFT arises from the connection between the external and 
KS potentials with the functional derivatives of the ki- 
netic and internal energy functionals. For CDFT these 
connections take the form (neglecting, for simplicity, 
again the spin degrees of freedom) 



ST, 



Sn{r) 



1 



A.(r)^ 



and 



Sip{r) c 



(17) 



(18) 



i.e., derivatives of the noninteracting kinetic energy 
determine the KS potentials Vs and As. Similarly, the 
derivatives of the internal energy functional defined 
as the ground-state expectation value of T and U, deter- 
mine the external potentials v{r) and A(r). Analogous 
equations hold also in DFT for superconductors. 

From Eqs. (|l^ and (|l|) we see that nonuniqueness of 
the Kohn-Sham potentials implies that the derivatives on 
the left-hand side do not exist on the space of all den- 
sities, for, if they existed, they would determine the po- 
tentials uniquely. Consequently, the functionals and 
F display multiple derivative discontinuities and must be 
redefined on equivalence classes of densities arising from 
the potentials modulo the nonunique pieces. The same 
applies to the xc functional E^c itself, since E^c is in gen- 
eral defined as the difference E^c = F — — Eh, where 
Eh stands for all Hartree-like terms included in the re- 
spective formulation of DFT. Common approximations 
to Ts and E^c do not display these derivative discontinu- 
ities. Judging from experience with similar discontinu- 
ities in ordinary DFT |16[| we expect this shortcoming to 
be most relevant for the calculation of energy gaps. 

The nonuniqueness problem discussed above (as well 
as the intimately related nondifferentiability problem) 
occurs, strictly speaking, only at zero temperature. At 
finite temperature one should work with a statistical en- 
semble, rather than with a ground-state, and then the 
uniqueness of the relation between density and potential 
is restored [0. However, the singularity at T=0 is an 
indicator that a real physical problem exists. Consider, 
as an illustration, the nonuniqueness of the potentials of 
SDFT, discussed in Refs. H and g. Due to the nondif- 
ferentiability of the xc functional, an infinitesimal change 
in the spin density (such as the change Sm caused by the 
flipping of a single electron in an extended half-metallic 
ferromagnet) may cause a finite (discontinuous) change 
in the xc potential. None of the existing approximations 
is able to reproduce such a discontinuity. Going to finite 
but small temperatures simply replaces the discontinu- 
ity by a very rapid continuous change. To estimate the 
scale of this change we note that at T = the mag- 
netic field is only determined by the densities to within 



Eg I fJ-o , where Eg is the energy gap and fiQ the Bohr mag- 
neton. Multiplying this with the low-temperature spin 
susceptibility we find that the spin density changes by 
{SmEG/kBT)exp{-EG/kBT). As long as keT « Eg 
this is much less than the physically relevant change 6m, 
and therefore the functional remains effectively discon- 
tinuous in the low-temperature regime. 

In summary, we have shown that generalizations of 
DFT to current-carrying and to superconducting sys- 
tems suffer from the same nonuniqueness problem we 
earlier discussed for the case of spin-polarized systems. 
Although the details are interestingly different in each of 
these three cases, the physical connection of nonunique- 
ness with a vanishing response function, as well as the 
classification of nonuniqueness into systematic (arising 
from constants of motion) and accidental (arising from 
special features of the ground state), and the conse- 
quences for differentiability of the respective density 
functionals and E^c are the same in all three cases, 
and, we believe, also in any other generalization of DFT. 
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